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ABSTRACT 

In the giant impact theory for lunar origin, the Moon forms from material ejected by the impact 
into an Earth-orbiting disk. Here we report the initial results from a silicate melt-vapor equilibrium 
chemistry model for such impact-generated planetary debris disks. In order to simulate the chemical 
behavior of a two-phase (melt+vapor) disk, we calculate the temperature-dependent pressure and 
chemical composition of vapor in equilibrium with molten silicate from 2000 to 4000 K. We consider 
the elements O, Na, K, Fe, Si, Mg, Ca, Al, Ti, and Zn for a range of bulk silicate compositions 
(Earth, Moon, Mars, eucrite parent body, angrites, and ureilites). In general, the disk atmosphere 
is dominated by Na, Zn, and O2 at lower temperatures (< 3000 K) and SiO, O2, and O at higher 
temperatures. The high-temperature chemistry is consistent for any silicate melt composition, and we 
thus expect abundant SiO, O2, and O to be a common feature of hot, impact-generated debris disks. In 
addition, the saturated silicate vapor is highly oxidizing, with oxygen fugacity (/o 2 ) values (and hence 
H2O/H2 and CO2/CO ratios) several orders of magnitude higher than those in a solar-composition 
gas. High fo 2 values in the disk atmosphere are found for any silicate composition because oxygen is 
the most abundant element in rock. We thus expect high oxygen fugacity to be a ubiquitous feature 
of any silicate melt- vapor disk produced via collisions between rocky planets. 

Subject headings: Moon — planets and satellites: composition — planets and satellites: formation — 
planet-disk interactions 



> 
o 

On 
cn 

o 

CO 



1. INTRODUCTION 

In the giant-impact theory, the Moon formed as a re- 
sult of an oblique collision between a Mars-sized impactor 
and the early Earth during its final stages of accr etion 
(jHartmann fc Davis! fl975l : I Cameron fe Wardl 11970 ). In 
the aftermath of this event, the Moon forms from the 
circumterrestrial disk of planetary debris ejected by the 
impact. The lunar composition thus provides constraints 
on the initial composition and subsequent chemical evo- 
lution of the protolunar disk. Indeed, a longstanding 
strength of the giant impact scenario is its ability to ac- 
count for a relatively Fe-poor, volatile-depleted Moon as 
an expected outcome of a large impact event with the dif- 
ferentiated proto-Earth (e.g., Hartmann & Davis 19751 : 
iStevensonl Il987l iCanup fe AsphaudfeOQll lCanupll2008l ). 

Dynamical simulations of the giant impact have com- 
monly found that the impactor pr ovides most (60%- 
80%) of the disk material (e.g., ICanup fe Asphaugl 
12001b iCanupl 12004 I2008D . However, major element 
abundances in the bulk silicate Moon closely resem- 
ble those in Earth's mantle, albeit with slight enrich- 
ments in refractory elements and sli ght depletions in 
volatile elem ents (e.g., see F i gure H IRingwoodl 119791 : 
iDrakd [1981 iRingwood et al.l 119871 ). Moreover, the 



Earth and Moon possess appare ntly identical isotopic 
compositions of O, Cr, and Ti (Wiechert et al.l 120011 : 
iLugmair fe Shukolvukovll 19981 : IZhang et al.ll2012f h which 
may require extensive re-mixing and equilibration be- 
tween the early Earth and the impactor-dominated ma- 
terial of the disk in order to remove any Earth-Moon 



differences (jPahlevan fe Stevenson! [2Q07h . 

Recent impact models suggest that a disk with a com- 
position essentially equal to that of the post-impact 
planet can result for scenarios that include a rapidly 
rotating early Earth dCuk fe Stewart i2012f ) or a suffi- 
ciently large impactor (jCanupl l2012f ). In these scenar- 
ios, Earth-Moon differences could presumably arise from 
the s ubsequent chemical evolution of the protoluna r disk 
(e.g.. IPahlevan et al.ll20lTI : iDesch fe Taylorl 1201 1[ ). For 
example, the observed fractionation of Zn isotopes in lu- 
nar samples suggests volatile loss during lunar formation 
(jPaniello et al.ll2012D . 

In any case, the impact-generated debris is ex- 
pected to form a two-phase silicate disk consist- 
ing of a melt layer s u rrounded by a vapor atmo- 
sphere (IStevensonl H982_KomDson & Stevenson 119881 : 
IPahlevan fe S tevenson! 120071 : IWard 20 ill). However, the 
chemistry of the disk atmosphere, and its implications 
for the subsequent chemical evolution of the disk, re- 
mains poorly characterized, and only t he evolution of 
the F eO/MgO ratio has been modeled (Pahlevan et al 
l2011f) . Her e we thus extend the mo dels of iFeglev et al 
(|2012l ) and IFeglev fe Schaeferl (|2012[ ) in order to explore 
the melt- vapor equilibrium chemistry of the protolunar 
disk. In particular, we examine the chemical composi- 
tion of the disk atmosphere, vapor pressure behavior of 
the silicate melt, and the prevailing oxygen fugacity of 
the disk. 

This letter is organized as follows. We describe 
our thermochemical equilibrium calculations for differ- 
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Fig. 1. — Elemental abundances in the bulk silicate Moon 
(gray lines with open symbols) and bulk silicate Earth 
(BSE; black lines with filled symbols) for elements in- 
cluded in the code calculations (in order of atomic num- 
ber). Lunar abund ances are taken from multiple liter- 
ature sourc e s (e.g.. iMorgan et al.l [19781: iRingwoodl | 1979| : 
Tavlor 1982; Wanke fc Dreibus 1982; Ringwood et al. 1987; 
Buck fc Toksoz 1980; Delano 1985; Jones fc Delano 1989; 
O'Neill 119911 : iWarren 20051: lLodders fc Feglev gOllh : BSE 
abundances are taken from lPalme fc O 'Neill (20031). 

ent melt compositions in §2. In §3 we present our initial 
model results, including the chemical speciation of the 
saturated silicate vapor, the total vapor pressure over 
silicate melts, the oxygen fugacity, and H2O/H2 and 
CO2/CO ratios in the disk. We conclude in §4 with a 
brief summary of our results and general implications 
for impact-generated silicate debris disks. Although our 
models focus on conditions most relevant for lunar for- 
mation, our results are generally applicable for impact- 
generated disks following collisions between terrestrial 
planets. The major gases that form in such collisions may 
thus provide chemical signatures of hot, impact-formed 
disks in extrasolar planetary systems. 

2. COMPUTATIONAL METHODS 
2.1. Chemical Model 

We used the MAGMA code ([Feglev fe Cameron! fl987l : 
iSchaefer fe Feglevl 12004 ) to calculate the temperature- 
dependent chemical equilibrium composition of silicate 
melt and coexisting vapor. The MAGMA code is a 
robust and rapid mass-balance, mass action algorithm 
based upon the ideal mixing of complex components 
(IM CC) model developed b y Hastie et al. (see references 
in lFeglev fe Cameron! 1 19871 ). A complete description of 
the MAGMA code, its method of calculation, and ex- 
tensive validation against exp erimental data is g i ven b y 
iFedev fe Cameron! (fl987h and ISchaefer fe Feglevl (|200l . 
Zinc thermochemistry was recently added to the code cal- 
culations in order to explore the chemical behavior and 
potential volatile loss of Zn-b earing gases from th e pro- 
tolunar disk atmosphere (e.g.. iPaniello et aT1l2012f ). 

The MAGMA code was used to simultaneously solve 
for activity coefficients in the melt, partial pressures in 
the gas phase, and melt- vapor equilibria for the elements 
Si, Mg, Fe, Al, Ca, Na, K, Ti, Zn, and O, consider- 
ing 53 melt species and 32 gas species (see Appendix) 
from 1800-6000 K. Neither H nor C are presently con- 
sidered in the MAGMA code. However as noted by 
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TABLE 1 

Bulk Silicate Compositions used in Model Calculations 



oxide wt% 


BSE 


Moon 


Mars 


EPB 


Angrite 


Ureilite 


Si0 2 


45.40 


44.60 


45.39 


46.26 


40.32 


43.65 


MgO 


36.76 


35.10 


29.71 


31.55 


11.10 


39.98 


A1 2 3 


4.48 


3.90 


2.89 


3.27 


11.45 


0.37 


Ti0 2 


0.21 


0.17 


0.14 


0.16 


1.12 


0.10 


FeO 


8.10 


12.40 


17.21 


14.82 


20.10 


14.60 


CaO 


3.65 


3.30 


2.36 


2.58 


15.89 


1.23 


Na 2 


0.349 


0.050 


0.983 


0.110 


0.022 


0.065 


K 2 


0.031 


0.004 


0.111 


0.009 


0.007 


0.022 


ZnO 


6.7e-3 


2.0e-4 


7.6e-3 


4.0e-5 


5.0e-4 


0.027 


Total 


99.0 


99.5 


98.8 


98.8 


100.0 


100.0 


Mg mol% 


89.0 


83.3 


75.5 


79.1 


49.6 


83.0 



Note. - BSE: bulk sili cate Earth (|Palme & O'Neill [20031 ): 
Moon: bulk silicate Moon ([ O'Neill 1991); Mars: bulk silicate 
Mars (Lodders 2000: [Lodders fe Feglev 2011); EPB: Eucrite Par- 
ent Body, presumably 4 Vesta (see Lodders 200^); Angrite is 
the unweighted mean composition of D'Orbigny, Sahara 99555, 
Angra dos Reis, LEW 86010, an d LEW 87051 angrit es com - 
puted from [Wa rren fe Kalle mevnl (| 19901 ); IWarren et al.l (|1995l) : 
Mittlefehldt et al. (2002); Ureilite is the unweighted median com- 
position of the ALH 84136, ALHA 77257, Dinngo Pup Donga, 
Dyalpur, Havero, Goalpara, North Haig, Novo Urei, an d PCA 
82506 ureilit e s computed f rom Ringwood (1960); Vdovvkin ( 19701) : 
rWiikl {1972); Uarosewichl (fl990): Klock et al.l (119981 ). All abun- 
dances are renormalized to 100% for code calculations. 

(Fegl ev fc Sch aefer 2012), the high temperature solubili- 
ties of H and C in molten silicates are expected to be low 
based on molecular dynamics simulations of CO2 solu- 
bility and experimental measurements of water solubil- 
ity versus temperature. Lower temperatures were not 
considered because the MAGMA code is designed for 
comp letely molten systems ab ove the liquidus tempera- 
ture (Fe glev fc Cameronl [T987). Here we focus on the re- 
sults from 2000-4000 K, as this encompasses the expected 
range of temperatures in a two-phase protolunar disk 
from the photosphere to the midplane. These tempera- 
tures are maintained over the ~100 yr lifetime of the disk 
via a thermal feedback between viscous spreading and va- 
porization of the disk melt layer ([Thompson fc Stevenson! 
119881 : Ward 2012). Chemical equilibrium calculations are 
justified because equilibri um is reached fast e r tha n the 
disk cools or dissipates. Feg lev fc Schaeferl ((2012) cal- 
culated chemical lifetimes of 10 -2 to 10 -4 seconds for 
the oxidation of Si, Mg, or Fe atoms to their respective 
monoxide gases via reactions such as M + O2 — >• MO + O 
in a silicate vapor at 2000 K, demonstrating that chemi- 
cal equilibrium is quickly reached between the major Fe, 
Mg, O, and Si gases in the disk atmosphere. 

2.2. Melt Composition 

We considered several different melt compositions, in- 
cluding bulk silicate compositions for the Earth, Moon, 
and Mars. Dating by three relative chronometers ( 182 Hf- 
182 W, 26 Al- 26 Mg, and 53 Mn- 53 Cr) indicates that the 
parent body of the eucrite, howardite, and diogenite 
achondrites differe ntiated during; the first 1-10 Ma of so- 
lar system history dJacobsen et al. 2008; T rinquier et al.l 
120081 : iSchiller et all 1201 ll ). We thus also included the 
major groups of differentiated achondrite meteorites in 
our modeling by considering the average composition 
of the eucrite parent body (EPB; presumably 4 Vesta), 
angrites, and ureilites. The oxide abundances and data 
sources for each melt composition are listed in Table [U 




Fig. 2. — Chemical composition of the saturated silicate vapor as a function of temperature, in equilibrium with melt compo- 
sitions corresponding to a) bulk silicate Earth (BSE), b) bulk silicate Moon, c) bulk silicate Mars, d) the Eucrite parent body 
(EPB), e) angrites, and f) ureilites. Table [1] lists the oxide abundances for each case. 



To a rough approximation, the bulk silicate Earth and 
bulk silicate Moon can be assumed to represent the ini- 
tial and final composition, respectively, of the protolunar 
disk material. We studied a range of bulk lunar composi- 
tions from the literature (e.g., see Figure [1]), and adopted 
IQ'Neill (jl99lh values as our nominal bulk silicate Moon 
composition for the model results reported here. 

As we discuss below, differences in the model results 
are primarily caused by differing abundances of volatile 
elements. For example, the bulk silicate Earth has Na, 
K, and Zn abundances of 2590, 2 60, and 54 /ig/g, re- 
spectively (|Palme fc O'Neilll [2003). whereas the corre- 
sponding lunar abundances are generally a magn itude 
lower , at 260, 31, and 1.9 /ig/g, respectively ([O'Neilll 
119911 ). Model lunar compositions with intermediate 
abundances of volatiles yield chemical trends and vapor 
pressures intermedi ate to the "e nd-member" results for 
bulk silicate Moon (IQ'Neillll 19911) and bulk silicate Earth 
(jPalme fc Q T NeIIlll2QQ3h . Moreover, model compositions 
with relati vely lower Fe abundances and higher Mg num- 
bers (e.g.. IWarrenll2005| ) give similar results as our nom- 
inal lunar model because Fe and Mg have less influence 
on the overall chemical behavior of the disk atmosphere. 

3. RESULTS 

3.1. Chemical speciation of saturated silicate vapor 

The chemistry of the saturated silicate vapor (i.e., the 
disk atmosphere) is shown in Figure [2] for different melt 
compositions over the temperature range 2000-4000 K. 



A number of general trends are observed for all silicate 
melt compositions. At high temperatures (> 3000 K), 
the disk atmosphere is dominated by SiO, O2, and O, 
whereas more volatile species such as Na and Zn tend to 
dominate the silicate vapor at lower temperatures. The 
mole fraction abundance of O2 is relatively temperature- 
independent, and O2 comprises 20-30% of the silicate 
vapor for all melt compositions. 

The chemical behavior of a given element in the silicate 
vapor is a function of both its abundance in the melt and 
its saturation vapor pressure behavior over the melt. By 
definition, volatile elements are characterized by higher 
vapor pressures. Thus, at lower temperatures, refrac- 
tory elements preferentially remain in the melt whereas 
volatile elements more readily enter the vapor, so that 
Na- and Zn-bearing species dominate the disk atmo- 
sphere. As temperatures increase, however, more disk 
material is vaporized and even refractory elements (e.g., 
Ca, Ti, Al) enter the vapor in increasing amounts. 

At high temperatures, the relative elemental abun- 
dances in the system play a larger role in determining 
the vapor composition, so that silicon (the most abun- 
dant element in the system after oxygen) will dominate 
the vapor as SiO, followed by O2, O, and the moder- 
ately volatile Mg- and Fe-bearing gases, which typically 
comprise 1%-10% of the vapor (e.g., Fe-oxides in the 
melt dissociate into Fe, FeO, O, and O2 in the vapor; 
iKazenas fc Tagirovl [l995). The temperature of the tran- 
sition between a SiO-dominated regime and a Na- and/or 




Fig. 3. — Model results showing a) total vapor pressure of the saturated silicate vapor for differ ent melt compositions, and b) 
oxygen fu gacity as a function of temperature. Al so shown are r esults for pure molten SiCb (e.g., [Thompson fc Stevensoni fl988; 
Ward 2012), and fo 2 for a solar-composition gas (Lodders 2003) and common oxygen fugacity buffers. MH: magnetite-hematite, 
QFM: quartz-fayalite-magnetite, IW: iron-wiistite, QFI: quartz-fayalite-iron. 



Zn-dominated regime depends upon the relative abun- 
dances of the volatiles in the melt. For example, for 
bulk compositions with relatively higher abundances of 
volatiles (BSE, Mars, Ureilite) this transition occurs be- 
tween 3000-3500 K. For melt compositions with lower 
volatile abundances (Moon, EPB, Angrite), this transi- 
tion occurs at temperatures near 2500 K. Moreover, as we 
discuss below, these basic differences in volatile content 
and disk chemistry also strongly influence the total vapor 
pressure and oxygen fugacity of the disk atmosphere. 

3.2. Total vapor pressure of saturated silicate vapor 

The vapor pressure of the saturated silicate vapor 
in equilibrium with the melt is an important param- 
eter that describes the therm al str ucture of the disk. 
iThompson fc Stevenson] (|1988f ) and IWardl (|2012[ ) used 
the Clausius-Clapeyron relation to describe the pressure- 
temperature profile of a two-phase disk, 

P = P oe -(T./T) } (1) 

and adopt iKriegerl ([1965) data for vapor pressure over 
pure molten Si02, using values of P Q = 2.23 x 10 8 
bar and T Q = 6.03 x 1 Q 4 K in the above expression. 
iFeglev fc Schaeferl (|2Q12h re-examined the vapor pressure 
behavior of molten silica up to 6000 K by considering 
chemical equilibria between molten silica and all species 
in the saturated vapor using the MAGMA code, giving 
the expression: 

logPva P (Si0 2 ,liq) = 8.203 - 25898.9/T, (2) 

for P in ba rs and T in K (cf. Melosh 2007]). Ho wever, 
as noted by IFeglev fc Schaeferl (|2012[ ) , the lunar disk is 
not pure silica but instead a mixture of multiple oxide 
components. We therefore calculated the total satura- 
tion vapor pressure in equilibrium with multicomponent 
silicate melts. The results are shown in Figure [3^. 

The saturation vapor pressure over a multicomponent 
melt is the sum of the partial pressures of all gases in 
equilibrium with a melt containing multiple oxides, 

Pt = SPi, (3) 



where i is a gas species (see Appendix) with a partial 
pressure Pi that is a function of the melt- vapor equilibria. 
For example, the net thermochemical reaction for SiO 
vapor can be written as 

Si0 2 (liq) = SiO(g)+0.5O 2 (g). (4) 

The partial pressure of SiO in the vapor is given by 

Psio = ilpsio 2 /o 2 ' 5 (5) 

where is the temperature-dependent equilibrium con- 
stant, fo 2 is the oxygen fugacity, and asio 2 is the ther- 
modynamic activity of liquid Si02 in a multicompo- 
nent melt. For volatile-bearing silicate melts, the sat- 
urated vapor is dominated by Na- and Zn-bearing gases 
at lower temperatures and by SiO at higher tempera- 
tures. Melt compositions with relatively higher volatile 
element abundances (e.g., BSE, Mars, Ureilite) therefore 
tend to have relatively higher saturation vapor pressures 
(Pt) at lower temperatures, as is demonstrated in the 
vapor pressure curves in Figure [3fci. The vapor pressure 
(bar) as a function of temperature is approximated by 

\ogP vap (BSE) « -51.17- 1048.8/T+ 14.60 logT (6) 

for a BSE-composition melt (jPalme fc O'Neillll2003[ ). and 

logP„ ap (BSM) w -41.45 -4052.3/T+ 14.90 log T (7) 

for a BSM-composition melt (jO'Neilll fl99lh , over the 
range of temperatures relevant for the disk atmosphere 
(2000-4000 K). 

Our models do not currently include H2O, and a 
full treatment of water solubility in high-temperature 
(> 1800 K) silicate melts is beyond the scope of this 
letter. However, we are able to make a number of gen- 
eral predictions regarding the behavior of water in the 
disk. Because of the relatively low solubility of water 
in magma (especially for high-temperature, low-pressure 
post-impact conditions), most H2O in the disk will re- 
main in the gas phase. Thus, for even minor amounts of 
water in the disk, we expect H2O to be a major gas in 
the disk atmosphere (for similar reasons that Na and Zn 
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Fig. 4. — Model results showing a) H2O/H2 ratios, and b) CO2/CO ratios in sa turated silicate vapor as a function of temper- 
ature. The corresponding ratios for pure molten SiCb, a solar-composition gas (|Lodd ers 2003), and common oxygen fugacity 
buffers are shown for comparison. 



dominate at lower temperatures), lead ing to an increase 
in the total vapor pressure of the disk. iFeglev fc Schaeferl 
(2012) reached a similar conclusion for the post-impact 
silicate vapor atmosphere of the Earth. We therefore con- 
sider any vapor pressure estimates for anhydrous melts 
(whether pure molten Si02 or silicate magmas) to repre- 
sent lower limits for the thermal structure of disk atmo- 
spheres. 

3.3. Oxygen fugacity of saturated silicate vapor 

The equilibrium chemistry between the silicate magma 
and coexisting vapor determines the overall oxygen fu- 
gacity (/o 2 ) of the system, given by 

fo 2 = lo 2 Po 2 = jo 2 Xo 2 P T , (8) 

where Xo 2 is the mole fraction abundance (see Figure [2j), 
Po 2 i s the partial pressure, and 7o 2 is the fugacity coef- 
ficient of O2 in a gas of total pressure Pt (see Figure [3^l) . 
Over the range of temperatures and pressures considered 
here, 7o 2 ~ 1, so we take fo 2 ~ Po 2 ~ Xo 2 P T from 
our model results. Figure [3}} shows the oxygen fugacity 
as a function of temperature for different silicate melt 
compositions, including pure molten silica. Also shown 
for comparison are fo 2 values for some common oxygen 
fugacity miner al buffers and for a solar-composition gas 
([Lodd ers 2003). Differences between model compositions 
at lower temperatures are again caused by the relative 
abundances of volatiles (e.g., Na, Zn) in each melt, and 
the fo 2 values converge at high temperatures where SiO 
is the dominant species in the disk atmosphere. For each 
case the saturated silicate vapor is highly oxidizing, with 
fo 2 values several orders of magnitude higher than those 
expected in a solar-composition gas. In fact, we expect 
very high /o 2 values for any silicate melt composition 
because oxygen is the most abundant element in rock. 

Although H2O and H2 are not included in the magma 
models, the prevailing oxygen fugacity — which is con- 
trolled by silicate melt-vapor equilibrium chemistry - 
can be used to constrain H2O/H2O ratios in the disk. 
For the net thermochemical reaction, 

H2+O.5O2 =H 2 0, (9) 



the equilibrium constant expression can be rearranged 
to give the H2O/H2 ratio as a function of the oxygen fu- 
gacity (/o 2 ) an d the temperature-dependent equilibrium 
constant 

H 2 0/H 2 = K$° l (10) 

where fo 2 is given by the silicate melt- vapor equilibrium 
chemistry. 

In the same way, the prevailing oxygen fugacity of the 
saturated silicate vapor can also be used to determine 
the relative abundances of CO2 and CO via the reaction, 

CO + 0.5O2 = C0 2 , (11) 
where the CO2/CO ratio is given by 

C0 2 /CO = Egtf^. (12) 

The H2O/H2 and CO2/CO ratios as a function of tem- 
perature are shown in Figures |4^l and HJd, respectively. 
Because these ratios are determined by the oxygen fu- 
gacity in the system, they are again significantly higher 
for an impact-generated disk atmosphere than for a solar- 
composition gas. For the silicate melts, the H2O/H2 ra- 
tio ranges from ~ 20 (Moon, Angrite, EPB) to ~ 100 
(Ureilite, Earth, BSE) at 1800 K to a nearly constant 
ratio of ~ 3 at higher temperatures (> 3000 K) for all 
silicate compositions. Similar behavior is observed for 
CO2/CO ratios, with values from ~ 6 (Moon, Angrite, 
EPB) to - 30 (Ureilite, Mars, BSE) at 1800 K to a nearly 
constant ratio of ~ 0.3 at higher temperatures. For com- 
parison, H2O/H2 - 0. 0003 and CO2 /CO < 0.0001 in a 
solar composition gas (Lodders 2003j). 

4. SUMMARY AND CONCLUSIONS 

We used thermochemical equilibrium calculations to 
model the chemistry of the circumterrestrial melt- vapor 
disk resulting from the Moon- forming impact. Our re- 
sults show that the silicate vapor in the disk is domi- 
nantly composed of SiO, O2, O at higher temperatures 
and Na, Zn, and O2 at lower temperatures. Although 
similar chemical trends are observed for all model com- 
positions, the relative abundances of volatile elements in 
the melt (particularly Na and Zn) strongly influence the 
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chemical composition, vapor pressure, and oxygen fugac- 
ity of the disk atmosphere at lower temperatures. Al- 
though not yet incorporated into our models, we expect 
similar chemical behavior from other volatile elements 
(e.g., Rb, Cs, Pb, Sn, etc.) which will readily enter the 
disk atmosphere. The chemical evolution of the disk due 
to vo latile loss (e.g., for Zn as reported bv lPaniello et al.l 
l2Q12f ) will be explored elsewhere. 

For any silicate composition, the silicate vapor in equi- 
librium with the melt is significantly more oxidizing (by 
several orders of magnitude) than a solar-composition 
gas. High fo 2 values in turn indicate high CO2/CO 
ratios and high H2O/H2 ratios, suggesting that only 
minor amounts of H2 and H will be produced by H2O 



dissociation over the lifetime of the disk. Because oxy- 
gen is the most abundant element in rock, high oxygen 
fugacity is a key feature of the protolunar disk, or any 
impact-generated debris disk that forms via collisions 
between rocky planets. High abundances of SiO, O2, 
and O may thus serve as chemical fingerprints of hot, 
impact-formed disks in extrasolar planetary systems. 
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nar Origin and Evolution (CLOE) as part of the NASA 
Lunar Science Institute (NLSI). Work by BF was sup- 
ported by NASA cooperative agreement NNX09AG69A 
with the NASA Ames Research Center and by grants 
from the NSF Astronomy Program. 



APPENDIX 

SPECIES INCLUDED IN MAGMA CALCULATIONS 

Melt Species: Si0 2 , MgO, FeO, Fe 2 3 , Fe 3 4 , CaO, A1 2 3 , Ti0 2 , Na 2 0, K 2 0, ZnO, MgSi0 3 , Mg 2 Si0 4 , MgAl 2 4 , 
MgTi0 3 , MgTi 2 5 , Mg 2 Ti0 4 , Al 6 Si 2 0i 3 , CaAl 2 4 , CaAl 4 7 , Cai 2 Ali 4 33 , CaSi0 3 , CaAl 2 Si 2 8 , CaMgSi 2 6 , Ca 2 MgSi 2 7 , 
Ca 2 Al 2 Si0 7 , CaTi0 3 , Ca 2 Si0 4 , CaTiSi0 5 , FeTi0 3 , Fe 2 Si0 4 , FeAl 2 4 , CaAli 2 0i 9 , Mg 2 Al 4 Si 5 0i 8 , Na 2 Si0 3 , Na 2 Si 2 5 , 
NaAlSi0 4 , NaAlSi 3 8 , NaA10 2 , Na 2 Ti0 3 , NaAlSi 2 6 , K 2 Si0 3 , K 2 Si 2 5 , KAlSi0 4 , KAlSi 3 8 , KA10 2 , KAlSi 2 6 , K 2 Si 4 9 , 
KCaAlSi 2 7 , Zn 2 Si0 4 , ZnTi0 3 , Zn 2 Ti0 4 , ZnAl 2 4 

Gas Species: O, 2 , Mg, MgO, Si, SiO, Si0 2 , Fe, FeO, Al, AlO, A10 2 , A1 2 0, A1 2 2 , Ca, CaO, Na, Na 2 , NaO, Na 2 0, Na + , K, 
K 2 , KO, K 2 0, K+, Ti, TiO, Ti0 2 , Zn, ZnO, e" 
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